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ABSTRACT 

We have derived the uncertainties to be expected in the derivation of galaxy physical properties 
(star formation history, age, metallicity, reddening) when comparing broad-band photometry to the 
predictions of evolutionary synthesis models. We have obtained synthetic colors for a large sample 
(~9000) of artificial galaxies assuming different star formation histories, ages, metallicities, reddening 
values, and redshifts. The colors derived have been perturbed by adopting different observing errors, 
and compared back to the evolutionary synthesis models grouped in different sets. The comparison has 
been performed using a combination of Monte Carlo simulations, a Maximum Likelihood Estimator and 
Principal Component Analysis. After comparing the input and derived output values we have been 
able to compute the uncertainties and covariant degeneracies between the galaxy physical properties as 
function of (1) the set of observables available, (2) the observing errors, and (3) the galaxy properties 
themselves. In this work we have considered different sets of observables, some of them including the 
standard Johnson/Cousins (UBVRcIc) and Sloan Digital Sky Survey (SDSS) bands in the optical, 
the 2 Micron All Sky Survey (2MASS) bands in the near-infrared, and the Galaxy Evolution Explorer 
(GALEX) bands in the UV, at three different redshifts, 2=0.0, 0.7, and 1.4. This study is intended to 
represent a basic tool for the design of future projects on galaxy evolution, allowing an estimate of the 
optimal band-pass combinations and signal-to-noise ratios required for a given scientific objective. 

Subject headings: galaxies: photometry - galaxies: evolution - methods: numerical - methods: 
statistical 



1. INTRODUCTION 

The catalogs produced by wide-field and all-sky sur- 
veys currently under development (e.g., GALEX; Martin 
et al. 1997, SDSS; York et al. 2000, 2MASS; Skrutskie 
et al. 1997, DENIS; Epchtein et al. 1997) in combination 
with astronomical databases like the NASA/IPAC Extra- 
galactic Database (NED) are beginning to provide easy 
access to extensive cross-correlated UV, optical, and near- 
infrared (NIR) photometry for millions of galaxies. The 
comparison of this huge amount of photometric data with 
the predictions of state-of-the-art galaxy population syn- 
thesis models provides an opportunity to obtain a more 
complete picture of the evolution of galaxies. One obvious 
goal is to gain insight into the star formation history, as 
well as the chemical and dust content evolution of galaxies 
from the high and intermediate redshift Universe down to 
the present. 

However, the reliability and precision of the derived 
galaxy properties expected to be found from these studies 
will depend on many factors; among these are: (1) the 
number of bands and wavelength coverage of the available 
observations, (2) the measuring uncertainties, and (3) the 
degeneracies between the different galaxy properties given 
the available photometric bands. In order to discriminate 
between different scenarios of galaxy evolution, the com- 
parison of photometric data and evolutionary synthesis 
models should also include the quantification of the un- 
certainties and covariances between the galaxy properties 
derived. In addition, the use of different external inputs 



(evolutionary tracks, stellar atmospheres libraries, etc.) in 
the evolutionary synthesis models leads to various discrep- 
ancies in the output results which may also result in sig- 
nificant uncertainties in the properties derived. While, a 
detailed study of these latter effects is beyond the scope of 
this paper, a nice approach to this problem may be found 
in Chariot, Worthey & Bressan (1996; see also Bruzual 
2000, Ccrvino, Luridiana & Castander 2000, Cerviho et 
al. 2001). 

Some recent studies have incorporated the effects of the 
observing errors and flux calibration uncertainties in the 
determination of the properties of different galaxy samples 
(Gil de Paz et al. 2000a, GILOOa hereafter; Bell & de Jong 
2000; Brinchmann & Ellis 2000). However, this type of 
analysis has not yet been performed in a systematic way, 
covering a large range of galaxy properties and/or star for- 
mation scenarios. A relevant exception is the work of Ro- 
nen, Aragon-Salamanca, & Lahav (1999) on the Principal 
Component Analysis of synthetic galaxy spectra. Simi- 
larly, no significant effort has yet been devoted to deter- 
mining quantitatively the optimal set of observables and 
signal-to-noise ratio required to obtain reliably derived 
galaxy properties. Noteworthy exceptions are the works 
of Kodama, Bell & Bower (1999), Bolzonella, Miralles & 
Pello (2000), and Wolf, Meisenheimer & Roser (2001) with 
regard to the galaxy classification and redshift determina- 
tion in broad and medium-band surveys. 

In this paper we explore the uncertainties expected in 
the derived properties of galaxies obtained from the analy- 
sis of broad-band photometry of nearby, intermediate, and 
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high-redshift objects. Wc quantify these uncertainties us- 
ing a sample of artificially generated spectral energy distri- 
butions (converted to broadband colors) with known (in- 
put) physical properties, studying the dependence of the 
properties derived with the sets of bands available, the ob- 
serving errors, as well as the galaxy star formation history 
and redshift. This work has interesting predictive capabil- 
ities and is intended to help in the design of future projects 
on the study of galaxy evolution. It can be used in the 
optimization of observing programs, helping to select the 
best wavelength bands and signal-to-noise ratios required 
to derive precise galaxy properties. This contribution is 
mainly focused on the study of the integrated properties 
of galaxies, star clusters and Hn regions, although it is our 
intent to generalize the application to spatially resolved 
portions of galaxies as well. A similar approach has been 
followed by Chariot & Longhetti (2001) for the optimiza- 
tion of emission-line data in galaxy spectra, although no 
error analysis was carried out by those authors. 

In the Section 2 of this paper we briefly describe the 
synthetic galaxy sample. The procedure followed for de- 
termining the uncertainties and degeneracies between the 
physical properties in this sample is described in Section 3. 
Section 4 includes a detailed description of the results from 
this analysis. The main conclusions are given in Section 5. 
Finally, some future applications for this work are given 
in Section 6. 

2. THE SAMPLE 

We have generated a large number of synthetic galaxy 
colors parameterized by different star formation histories, 
metallicities, reddening, and redshifts. Because of the un- 
bounded number of possible combinations we have made 
some choices and simplifications. 

With regard to their star formation history we have con- 
sidered galaxies with exponential star formation having 
timcscales (r) between 0.2 Gyr and 6Gyr. Although this 
scenario is a rough approximation to the actual star for- 
mation histories it constitutes the most widely accepted 
parameterization of the star formation history in individ- 
ual galaxies (see e.g. Brinchmann & Ellis 2000). 

Thus, the properties to be determined are the timescale 
for the galaxy formation, the age of the galaxy, the stel- 
lar metallicity, and the reddening. A total of 9000 syn- 
thetic galaxies were generated, a third of them at redshift 
2=0.0, another third at z=0.7, and the remainder at red- 
shift 0=1.4. These redshift values were chosen to cover 
the epoch where most of the evolution of the star forma- 
tion activity in the Universe has apparently taken place 
(Gallego et al. 1995; Madau et al. 1996; Connolly et al. 
1997; Madau, Dickinson & Pozzeti 1998). Since the mod- 
els available only provide discrete values in metallicity we 
chose to assign to each galaxy the nearest metallicity value 
of those given by the corresponding model. The ranges of 
physical properties covered by the galaxies in the sample 
are shown in Table 1. 

Although a complete study about the effects of the 
model uncertainties on the galaxy properties derived is 
beyond the scope of this paper, the impact of the errors 
in the broad-band colors due to stellar evolution prescrip- 
tion and spectral calibration uncertainties (Chariot et al. 
1996; Yi, Demarque & Oemler 1997) will be taken into 



account in our further analysis. In addition, in order to 
illustrate the effect of these uncertainties we have consid- 
ered two different sets of evolutionary synthesis models. 
The synthetic galaxies were generated using the predic- 
tions for the stellar continuum given by the GISSEL99 
models (Bruzual & Chariot, in preparation) while the 
best-fitting set of properties has been derived using both 
the GISSEL99 and the PEGASE models (Version 2.0; see 
Fioc & Rocca-Volmerange 1997). Although the theoretical 
isochrones in both models come mainly from the Padova 
group (Bressan et al. 1993), there are fundamental differ- 
ences in post-main-sequence evolutionary stages between 
the two set of models, including in the early, thermally 
pulsating, and post-AGB phases. Noteworthy the differ- 
ences in these evolutionary stages are responsible for the 
most serious uncertainties in the stellar population mod- 
eling (see Chariot et al. 1996). In Figure 1 we show the 
predictions of these two models for Simple Stellar Popu- 
lation (SSP) and continuous star formation galaxies with 
Solar metallicity. For the sake of comparison we have also 
included the predictions of the Starburst 99 models (Lei- 
therer et al. 1999) at ages younger than 0.1 Gyr. 

3. ANALYSIS 

3.1. Galaxy Colors 

Once the sample was generated we determined the lumi- 
nosity (per Solar mass) in the different bands and the col- 
ors for each individual galaxy in the sample using the pre- 
dictions of the corresponding evolutionary synthesis mod- 
els. We adopted for all the models the same Salpeter IMF 
with Mi ow =O.lM and M up =100M Q . 

The bands considered in this work include the John- 
son/Cousins (U BV Rc Ic) an d Sloan Digital Sky Survey 
u'g'r'i'z' (Fukugita et al. 1996) optical bands, the JHK 
near-infrared bands, and the GALEX near-ultraviolet 
(NUV, 1800-3000A) and far-ultraviolet bands (FUV, 1350- 
1800 A; Doliber et al. 2000). Because of the small differ- 
ences expected between the standard isT-band and the K s 
and K' bands, we decided to include only the standard 
K in our realizations. In this sense, although the 2MASS 
survey has been carried out using -Kg-band imaging data, 
we will refer hereafter to the 2MASS data set when dealing 
with the standard JHK bands. 

The colors obtained for each stellar population were 
then reddened (and the apparent mass-to-light ratios in- 
creased) using the corresponding E(B — V) values and 
adopting the parameterization of the Galactic extinction 
law of Cardelli, Clayton & Mathis (1989) for a total-to- 
selective extinction ratio i2y=3.1. However, since this pa- 
rameterization is only valid down to 1000A, and given that 
some of the bands selected for this work cover regions of 
galaxy spectra well below the Lyman limit (for z=0.7 and 
0=1.4), we extended the extinction law to the far-UV using 
the A\/Ay ratios given by Mathis (1990; see also Martin 
& Rouleau 1989). The A\/Ay mean values for each band 
(its respective redshiftcd rest-frame wavelength) were then 
computed convolving the filters response functions with 
the adopted extinction curve. 

Using the number of ionizing Lyman photons predicted 
by the evolutionary synthesis models we also computed 
the contribution of the nebular continuum and most in- 
tense gas emission-lines to all the bands considered. We 
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assumed that 85% of the photons with A<912 A effectively 
ionize the surronding gas, but a 15 per cent fraction would 
be observed at the far-ultraviolet or absorbed by dust (Lei- 
therer et al. 1995; Dove, Shull & Ferrara 2000; GILOOa). 

Most of the 85% of the far-UV photons absorbed by 
the surronding gas in the Hll region is reprocessed as 
nebular continuum and emission lines in the optical and 
NIR. However, a very small fraction of these photons can 
be re-emitted as free-free radiaton in the far-UV. In our 
case, we assumed that these secondary far-UV photons are 
not absorbed by neutral gas again, but rather being ab- 
sorbed by dust or escaping from the galaxy. In order to 
determine the nebular continuum contribution to all the 
bands we have used the emission and recombination coef- 
ficients given by Ferland (1980) for the near-UV, optical, 
and NIR for T o =10 4 K. For the far-UV free-free radiation 
we have assumed a constant gaunt factor g^ — 1.1 (see 
Karzas & Latter 1961) in the range 500-912 A for a gas 
with 7W0 4 K. 

With regard to the gas emission-lines we have 
assumed the relation between the number of Ly- 
man photons and Ha luminosity given by Brocklc- 
hurst (1971) and the theoretical hydrogen line-ratios 
expected for a low density gas (ra c =10 2 cm~ 3 ) with 
T C =10 4 K in Case B recombination (Osterbrock 1989). 
We considered the contribution of the most intense 
forbidden lines ([On]AA3726,3729A, [Om] AA4959,5007A, 
[Nn]AA6548,6583A, [Sn]AA6717,673lA) adopting the mean 
line ratios measured by Gallego et al. (1996) for the Uni- 
versidad Complutense de Madrid (UCM) sample of lo- 
cal star-forming galaxies (Zamorano et al. 1994, 1996). 
Both the nebular continuum and the emission-line lumi- 
nosities were corrected for extinction assuming the rela- 
tion given by Calzetti, Kinney & Storchi-Bergmann (1996, 
see also Calzetti 1997; Storchi-Bergmann, Calzetti & Kin- 
ney 1994) between the gas and the stellar continuum 
reddening associated with the young stellar population: 
E(B - V) stol iar=0.44x£(B - U) gas . 

Following this procedure we obtain all the input infor- 
mation concerning the actual properties, colors, and mass- 
to-light ratios for the galaxies in the sample. In the next 
section we discuss our recovery of these input properties 
starting from the observed galaxy colors, and their corre- 
sponding measuring errors. 

3.2. Sets of Colors 

Grouping the colors deduced for these galaxies in differ- 
ent sets and comparing them with the evolutionary synthe- 
sis models allows us to explore systematics and determine 
the set of observables that result in minimum differences 
between the actual (input) galaxy properties and the de- 
rived (output) properties. This comparison is done using 
a combination of Monte Carlo simulations, a maximum 
likelihood estimator, and a Principal Component Analysis 
algorithm. 

The number of different combinations of colors that 
could be constructed considering a total of 10 potential 
bands from the UV to the near- infrared is X^r=2 r\(w-r)\ = 
1013. For this study we have selected only the 10 sets 
shown in Table 2. The detailed comparison of the results 
obtained for all these sets provides enough information 
about the relevance of the different bands, observing er- 



rors, etc., for a precise determination of the galaxy prop- 
erties. 

3.3. Comparison Procedure 

Once the colors of the galaxy sample had been obtained 
and grouped in sets we then perturbed the "observed" 
magnitudes by applying random observing errors. In or- 
der to simplify the problem we studied three cases, cor- 
responding to three different 1-sigma errors in the colors, 
0.03, 0.07, and 0.10 mag, and consider only the case where 
these errors are the same in all the colors. 

In order to compute the effects of these observing errors 
in the galaxy properties to be derived we used a Gaussian 
distribution of errors for all bands generated using a Monte 
Carlo simulation method. The colors derived for each test- 
particle were then compared with the evolutionary synthe- 
sis models using a maximum likelihood estimator, L. The 
expression for this estimator (see e.g. Abraham et al. 1999) 
is 

where C n are the colors derived, c„ are those predicted 
by the evolutionary synthesis models and N is the total 
number of colors available within each set. Because the 
same level of error was assumed for all the bands, the 
maximization of this expression is equivalent to comput- 
ing the minimum \ 2 ■ Therefore, in this case, we could 
estimate the confidence levels in the galaxy properties via 
the Avni's approximation (Avni 1976) as has been done 
by Bolzonella et al. (2000) instead of using numerical sim- 
ulations. However, the considerations that lead to this 
estimation procedure apply only asymptotically, being ap- 
plicable when the x 2 estimator covariance matrix can be 
replaced by its linear approximation in the vicinity of the 
best-fitting set of parameters. Although these conditions 
could be fulfilled in our case (Perez-Gonzalez et al. 2001, 
in preparation) , we have decided to use numerical simula- 
tions in order to be able to derive the degeneracies between 
the different galaxy properties. 

The ranges in the evolutionary synthesis models param- 
eters where the data-model comparison was performed are 
shown in Table 1. In order to avoid introducing a con- 
straint bias in the derived properties, the ranges for these 
comparison were chosen to be significantly wider than 
those where the galaxy sample was generated. 

Once the expression C is maximized for a significant 
number of Monte Carlo test-particles (we used a total of 
200) we obtained the distribution of physical properties 
associated with the probability distribution of the galaxy 
colors. In Figures 2a, b & c we show the results obtained 
for a nearby galaxy with a exponential star formation his- 
tory with t=4.5 Gyr timescale, Solar mctallicity, an age of 
5.0 Gyr, and E(B — V) of 0.08 mag, for observing errors in 
the colors of 0.10, 0.07, and 0.03 mag, respectively. In this 
case the set of colors used was that including the GALEX, 
SDSS, and 2MASS bands (see Table 2). 

These figures illustrate the strong degeneracy between 
the different galaxy properties even for relatively small 
measurement errors. We have therefore performed a quan- 
titative analysis of these degeneracies using a Principal 
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Component Analysis (PCA hereafter) on the space of 
galaxy properties by solving the eigenvector equation on 
the test-particle correlation matrix of each galaxy in the 
sample (see Morrison 1976). This analysis gives the direc- 
tion in the space of galaxy properties along which the 200 
solutions obtained for each individual galaxy are mainly 
oriented, constituting the best estimator of the degener- 
acy between these properties. 

Summarizing, this procedure provides us with (1) the 
mean derived properties, (2) the 1-sigma errors, (3) the 
orientation in the space of solutions of the Principal Com- 
ponent (PCAl hereafter), and (4) the input (i.e. actual) 
properties of all galaxies in the sample. 



4. RESULTS 

Once these quantities had been derived we computed 
the mean differences between the output and input val- 
ues along with the mean 1-sigma spread, at fixed inter- 
vals in the input properties. The bins used were 0.5 Gyr, 
0.025 dex, and 0.025 mag in the formation timescale, age, 
and reddening, respectively. Mean differences and 1-sigma 
values in the stellar mctallicity were computed for each of 
the input values considered. 

In Figures 3a & b we show the results obtained before 
and after computing the mean differences and 1-sigma er- 
rors for a subsample of 500 nearby galaxies with errors in 
the colors of 0.07 mag and U+BVRI+JHK data available. 
Due to the relevance of the i-T-band luminosities in order 
to derive stellar masses in nearby (Aragon-Salamanca et 
al. 1993; GILOOa) and intermediate-rcdshift galaxies (see 
Brinchmann & Ellis 2000 and references therein) , the mean 
differences between the derived and input iVband mass- 
to-light ratios were also computed. Mean uncertainties for 
all the sets of observables, redshifts, and observing errors 
considered are summarized in Table 3. 

In addition, we studied the degeneracies between the 
galaxy physical properties analyzing the distribution of 
the unitary PCAl vector components. In Figure 4a 
we show the frequency histograms obtained for the 
sample of high-rcdshift galaxies assuming an error in 
the colors of 0.07 mag and the SDSS, SDSS+2MASS, 
and GALEX+SDSS+2MASS sets. Note that the 
PCAl vector points toward the direction where the 
largest fraction of the galaxy properties' variance oc- 
curs. In this sense, a PCAl vector with compo- 
nents (ui g t ,u E (B-v),uiog z/z Q ,u T )=(+0. 707,-0. 707,0,0), 
say, implies the existence of a degeneracy between 
age and reddening in the sense that younger, ob- 
scured stellar populations have colors that are indis- 
tinguishable from older but less extincted ones. In 
this case, no age-metallicity or age-timescale degenera- 
cies would be present. However, the behavior described 
above could also result in a PCAl vector with compo- 
nents (ui g t,u E (B-v),uio g z/z ? ,u T )=(-0. 707, +0.707, 0,0). 
That, however, would appear in a different position at the 
frequency histogram shown in Figure 4a. Thus, in order 
to reduce this sign ambiguity when interpreting our results 
we forced the ui og t component to be positive, changing the 
sign of all the vector components if ui og t was negative. Fi- 
nally, in order to quantitatively determine the dominant 
degeneracy for each individual galaxy in our sample we 



have defined the angle 9ij like that satisfying 

cos 9 id = SIGN (^j y/v* + u) (2) 

where u, and Uj are the i and j components of the PCAl 
vector. The angle 6ij simultaneously provides a measure 
of the angle between the PCAl vector and the plane of 
physical properties i,j and the sign of the degeneracy be- 
tween the i and j properties. Thus, if | cos#j ; j|~l the de- 
generacy between the i and j properties would dominate 
the total degeneracy. Moreover, if cos9ij>0 an increase 
in both the i and j properties would lead to similar obser- 
vational properties, while if cos#jj<0 the value of one of 
the properties should decrease. In Figure 4b we show the 
distribution of cos 6ij as function of the age for the high- 
redshift sample assuming an error in the colors of 0.07 mag 
and the GALEX+SDSS+2MASS set available. 

Along this section we will describe the results obtained 
from the analysis of the distributions shown in Figures 3 
& 4 for the different redshifts, observing errors, and band- 
pass combinations considered. 

4.1. Nearby galaxies 

4.1.1. Formation Timescale 

With regard to the formation timescale in nearby galax- 
ies, Figure 5a indicates that, even for relatively small ob- 
serving errors, its uncertainty is very high (see also Ta- 
ble 3) and shows a strong dependence with the value of 
the formation timescale itself. The larger uncertainty in 
the formation timescale for larger values of this quantity 
is mainly due to the small sensitivity of the optical-NIR 
colors of stellar populations with ages t«r to changes 
in its formation timescale. The use of [/-band data sig- 
nificantly reduces this uncertainty, probably due to the 
high sensitivity of this band to the presence of recent star 
formation that allows to rule out instantaneous-burst so- 
lutions when recent star formation associated with larger 
t values has effectively taken place. The use of NIR data, 
however, does not provide relevant information about the 
formation timescale of the stellar population. Moreover, 
the reduction achieved in the uncertainties of the differ- 
ent galaxy properties by using JHK data compared with 
those obtained using exclusively iVband data is very small 
(see Table 3). As we will show in Sections 4.2 & 4.3 this is 
not the case for the intermediate and high-redshift galax- 
ies, where these bands now cover the redshifted optical 
spectrum. Finally, in the same way that the [/-band, the 
use of UV data provides an additional reduction in the 
formation timescale. As it is clearly seen in Figure 2, for a 
particular galaxy the formation timescale is mainly degen- 
erate with the age of the stellar population, in the sense 
that, within the observing errors assumed, an increase in 
the formation timescale accompanied by an increase in the 
age can result in similar UV-optical-NIR colors. Due to 
this age-timescale degeneracy part of the reduction in the 
timescale uncertainty obtained by the use of UV data can 
be explained by the significant reduction in the age uncer- 
tainty achieved by including UV data (see below). 

4.1.2. Age 

With respect to the age determination in nearby galaxies 
(z=0), the Figure 5a also shows that a significant reduc- 
tion in the age uncertainty is achieved by including NIR 
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data. It is important to note that the use of additional 
NIR data result in the same dominant degeneracy that 
if only optical data are used (see below), but the range 
of physical properties where this degeneracy takes place is 
significantly smaller. The most significant improvement in 
the age determination, however, is obtained when UV data 
are available. This is mainly due to the emission arising 
from post-AGB stars in low-metallicity populations and at 
ages younger than 10 Gyr and to the "UV-upturn" in high- 
metallicity evolved (i>10Gyr) stellar populations that re- 
sult in highly peculiar UV-optical colors. However, the 
uncertainty in the modeling of post-AGB stars (Chariot et 
a. 1996) and the low-mass core helium-burning Horizontal 
Branch (HB hereafter) and evolved HB stars that lead to 
the "UV-upturn" (Yi et al. 1997), introduce additional er- 
rors in the UV-optical colors during the data-models com- 
parison. Chariot el al. (1996) estimated using two differ- 
ent theoretical prescriptions that the uncertainty only in 
the post-AGB phase modeling could result in differences 
of about 1 mag in the UV-optical colors of a several-Gyr- 
old stellar population. Moreover, although we assumed 
the same observing errors for all the bands, the faint UV 
emission of evolved stellar populations is expected to re- 
sult in very large observing errors in the UV-optical colors. 
Therefore, while the stellar evolution of these stars is not 
well understood the age determination in old stellar pop- 
ulations should not rely on the use of UV data. 

Along with the formation timescale, the age of the stel- 
lar population in nearby galaxies is mainly degenerate with 
the dust extinction, in the sense that older stellar popula- 
tions with low dust content have similar colors to highly- 
extincted, younger stellar populations. Although in the 
case of very old stellar populations the age-extinction de- 
generacy also competes with the age-metallicity degener- 
acy (see Worthey 1994), the age-extinction degeneracy is 
still dominant in this range for all the band-pass combina- 
tions and observing errors considered in this work. It could 
be argued that the strong discretization of the metallicity 
in our models could be responsible for the relatively weak 
age-metallicity degeneration derived. However, the fact 
that this behavior is observed even for the largest observ- 
ing errors considered indicates that it is real and a natural 
consequence of the use of broad-band data. In this sense, 
the combination of broad-band with narrow- or medium- 
band data or spectroscopic indexes would break the age- 
extinction degeneracy, making of the age-metallicity the 
dominant degeneracy (see Worthey 1994). 

4.1.3. Dust Extinction 

The dust extinction is derived with a very high accuracy 
(E{B — V)=0.04-0.20mag) even for large observing errors 
and relatively low number of observables (see Table 3). 
In the case of the nearby galaxies, the uncertainty in the 
dust extinction does not depend on the value of the dust 
extinction itself and is mainly degenerate with the age of 
the stellar population (see above) with some contribution 
from the extinction-metallicity degeneracy. In combina- 
tion with BVRI optical data either the use of UV, U, 
or NIR data provide a significant reduction in the dust- 
extinction uncertainty. In order to better derive the dust 
extinction the use of a wider wavelength baseline in wave- 
length (e.g. using UVIJK) is more effective than fully cov- 



ering the optical range (UBVRI) . This is mainly due to the 
reduction in the metallicity uncertainty by the use of NIR 
data (see below) that leads, via the extinction-metallicity 
degeneracy, to a reduction in the dust-extinction uncer- 
tainty. Again, the use of JHK data instead of only K- 
band data do not lead to a significant reduction in the 
dust-extinction uncertainties. 

4.1.4. Metallicity 

With regard to the metallicity of the stellar popula- 
tion the uncertainties derived are strongly dependent on 
the band-pass combination available and the value of the 
metallicity itself. In particular, the uncertainties derived 
are smaller as the metallicity becomes higher (see Fig- 
ure 5a). Within the age range considered, the main con- 
tributors to the optical and NIR emission of SSP galaxy 
are the main-sequence and RGB stars. However, for a 
more constant star formation, a significant contribution 
from core-Helium-burning stars is expected (see Chariot 
& Bruzual 1991). In order to determine the source of 
the metallicity dependence of these uncertainties we have 
produced the same diagrams shown in Figure 5a but re- 
stricted to formation timescales shorter than 50 Myr. The 
analysis of this diagram shows no dependence of the un- 
certainties with metallicity, which implies that the source 
of the dependence was the distinct photometric evolution 
of high-metallicity core Helium burning stars (Mowlavi et 
al. 1998). It is worth noting, however, that at very high 
metallicities the uncertainties in the modeling of the stel- 
lar populations are themselves very large because of the 
lack of very metal rich stars of any age in the Solar neigh- 
borhood that could be used as spectral calibrators (see 
Chariot et al. 1996). 

The most significant reduction in the mean metallicity 
uncertainty is achieved when NIR data are used in com- 
bination with optical data (see Figure 5a). Although the 
uncertainties in the model predictions for the thermally 
pulsating AGB (TP-AGB hereafter) and the upper RGB 
can result in differences in the (V — K) color predicted by 
different models of -0.10-0. 15mag (Chariot et al. 1996), 
the improvement in the metallicity determination by the 
use of NIR data is still relevant. In this sense, in Ta- 
ble 3 we show that the mean metallicity uncertainty for 
the U+BVRI set is 0.32 dcx assuming an observing error 
of 0.03 mag, while the uncertainty for the U+BVRI+K set 
assuming an observing error of 0.10 mag is only 0.26 dex. 

4.1.5. Stellar Mass 

As input for the iVband mass-to-light ratio of the stellar 
populations we have adopted M#-.0=3.33 (Worthey 1994). 
It should be noticed that along with the errors in the stel- 
lar mass-to-light ratios derived here the misunderstanding 
of the actual IMF introduce an additional, systematic un- 
certainty, which, in fact, constitutes the most important 
source of error in the determination of the galaxy stellar 
mass (Bell & de Jong 2001). In addition, the poor con- 
straints on the theoretical isochrones of upper-RGB stars 
and AGB stars can result in a 20 per cent uncertainty in 
the -ftT-band mass-to-light ratio (Chariot et al. 1996). In 
Figures 6a & 6b we show the uncertainties expected in 
the if-band mass-to-light ratio from different sets of ob- 
servables that include if -band data. These uncertainties 
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show a strong dependence with the galaxy age and forma- 
tion timescale in the sense that larger uncertainties are ex- 
pected at lower values of the formation timescale and older 
ages. Figure 5a shows that the value of the age uncertainty 
(in log t scale) is almost independent of the age itself. In 
addition, Figure 1 indicates that the rate of change in the 
iVband mass-to-light ratio (with log t) is higher when the 
stellar population becomes older, specially for very low val- 
ues of the formation timescale. Therefore, for a constant 
uncertainty in log t, an increase in the uncertainty of the 
mass-to-light ratio at very old ages is expected. 

Figure 6a also shows that the mass-to-light ratio deter- 
mination is biased toward lower values. This bias, which 
is particularly important at old ages, is probably due to 
the upper limit of 15 Gyr in age imposed during the data- 
models comparison (see Table 1), although other contrib- 
utors can not be ruled out (see Sections 4.2.2 & 4.3.2). As 
it is clearly seen in Figure 6b, the use of UV data allows to 
reduce both the uncertainty and bias in the mass-to-light 
determination. This reduction is directly related with the 
reduction in the age uncertainty described above. How- 
ever, as we already commented, the use of UV data for the 
study of stellar populations with ages older than several 
Gyr can lead to wrong conclusions because of the uncertain 
modeling of the post-AGB phase and the "UV-upturn" . 

The behavior described above for the timescale, age, 
dust extinction, mctallicity, and stellar mass is identical 
for any observing error but with larger mean uncertainties 
for larger observing errors. The reader is referred to the 
Table 3 for the dependence of the mean uncertainties in 
the different galaxy properties derived with the observing 
errors. 

4.2. Intermediate-redshift galaxies 
4.2.1. Formation Timescale 

In Figure 7 we show the uncertainties derived for the 
properties of intermediate-redshift galaxies (z=0.7). With 
regard to the formation timescale the uncertainties are 
very large (2-3 Gyr), even larger than for the nearby galaxy 
sample. As we commented in Section 4.1, the optical- NIR 
colors are quite insensitive to changes in the formation 
timescale with t«T. Therefore, since we are assuming 
that these galaxies are statistically younger than the those 
observed in our Local Universe (see Table 1) and the range 
in formation timescale is obviously the same, the uncer- 
tainty in the formation timescale is necessarily higher. For 
the same reason the uncertainty at very low timescale val- 
ues is much lower than at high timescale values. 

The upper panel of Figure 7a also suggests a significant 
bias in the timescale determination toward lower values of 
this property. This bias is also the consequence of the small 
changes in the optical-NIR colors of these galaxies with 
the timescale when the age is younger than the timescale 
value. In this case, the higher rate of change in the colors 
toward lower formation timescales systematically leads to 
lower values in order to reproduce the probability distri- 
bution associated with the observing errors. It is worth 
noting that, because of the reduction of this bias, the use 
of a larger number of bands may result in some cases in 
a higher timescale uncertainty (see Table 3 for the results 
on the U+BVRI+K and U+BVRI+JHK sets). Like in 
the nearby galaxies case, the dominant degeneracy involv- 



ing the formation timescale is the age-timescale degener- 
acy, in the sense that older galaxies with high formation 
timescales have similar colors that younger galaxies with 
a more instantaneous star formation. This is true for any 
band-pass combination considered. With regard to the op- 
timal set of observables, Table 3 demonstrates that for the 
same number of bands the use of wider wavelength base- 
lines results in lower uncertainties. In particular, the use 
of the UVIJK set reduces the timescale, age, and metal- 
licity uncertainties inherent to the UBVRI set providing 
also a much lower dust-extinction uncertainty than the 
BVRI+K set. On the other hand, the SDSS+2MASS and 
GALEX+SDSS+2MASS sets result in very similar uncer- 
tainties (see Table 3), which implies that the optical and 
NIR bands provide most of the information available in 
the UV and in the blue part of the optical spectrum about 
the galaxy age, star formation history, and metallicity. 

4.2.2. Age 

With respect to age of the intermediate-redshift galaxies 
the uncertainties derived are significantly smaller than in 
the nearby-galaxies case. This is mainly due to the higher 
rate of change in the rest-frame optical colors within the 
age range assumed for these galaxies compared with that 
assumed for the nearby galaxies (see Table 1 and Figure 1). 
In addition, the fact that the iVband now corresponds to 
the rest-frame J-band emission implies that the effect of 
the uncertainties in the model predictions associated with 
the upper RGB and AGB evolutionary stages is less im- 
portant (see Section 4.4.2). On the other hand, the use 
of [/-band data for determining ages older than 1 Gyr at 
these redshifts is strongly limited by the uncertainty in 
the modeling of the rest-frame UV emission from post- 
AGB stars (Chariot et al. 1996; see Section 4.1). How- 
ever, the most significant decrease in the age uncertainty 
is achieved when NIR data are used, specially if data in 
all the bands (JHK) are available. This is probably due 
to the fact that the JHK set provides information simul- 
taneously about the presence of AGB stars (via the rest- 
frame z' and J bands) and main-sequence stars (via the 
rest- frame i?-band). 

Figure 5a also shows the existence of a significant bias 
toward younger ages for the BVRI and UBVRI sets. In 
this case the presence of this bias is due (1) to the existence 
of a formation timescale bias and a strong age-timescale 
degeneracy and (2) to the fact that the optical colors of the 
stellar populations change more slowly as the population 
becomes older. In the latter case, in order to reproduce 
the distribution of optical colors associated with the ob- 
serving errors, the best-fitting solution should be found at 
younger ages, where the intrinsic dispersion of the model 
colors is larger. As we show below a bias in age also results 
in a bias in the galaxy if -band mass-to-light ratio. The 
use of wider wavelength baselines allows to significantly 
reduce this bias. In particular, the use of the UVIJK leads 
to a less severe bias and lower age uncertainties than the 
U+BVRI and the U+BVRI+K sets. Within the age un- 
certainty interval the degeneracy is mainly dominated by 
the age-timescale degeneracy with some contribution from 
the age-extinction degeneracy in those band-pass combi- 
nations that do not include UV or [/-band data. 

4.2.3. Dust Extinction 
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The dust extinction in the sample of intermcdiate- 
rcdshift galaxies is derived with high accuracy specially 
when [/-band data are available (see Figure 7a). In 
this case, the uncertainties expected in E(B — V) are in 
any case smaller than 0.10 mag for observing uncertain- 
ties as high as AC„=0.10mag. The significant reduction 
achieved, if we compare these results with those derived 
for the nearby galaxies, is due to the very high sensitiv- 
ity of the redshifted UV emission to the presence of small 
amounts of dust. In those band-pass combinations not 
including [/-band data we notice a clear dependence of 
the dust-extinction uncertainty with the value of the ex- 
tinction itself, with larger uncertainties at larger values of 
the extinction (see Figure 7a). The analysis of the PCAl 
components also indicates that at dust-extinction values 
higher than E{B — V)>0.5mag the age-extinction degen- 
eracy becomes very important. This implies that in highly 
extincted intermediate-redshift galaxies a small increase in 
the amount of dust can lead to the same optical colors 
(specially if U -band data are not used) that a comparable 
decrease in the age of the stellar population would pro- 
duce. 

4.2.4. Metallicity 

With regard to the metallicity uncertainty, Figure 7a 
shows that the uncertainty decreases with the value of the 
metallicity itself. The reduction is particularly important 
when NIR data are available. The use of the three JHK 
NIR bands reduces this uncertainty over the whole range 
of metallicities. In this sense, the use of the UVIJK set re- 
sults in lower metallicity uncertainties than the U+BVRI 
and the U+BVRI+K sets (see Table 3). It is important 
to keep in mind that the JHK filters now cover the rest- 
frame R, z' , and J bands. In the age range considered 
the main contribution to the rest-frame optical emission 
comes from main-sequence stars. On the other hand, 
the rest-frame NIR emission, along with main-sequence 
stars, shows an important contribution from AGB and 
core-Helium-burning stars (see Chariot & Bruzual 1991). 
The role played by AGB stars is more relevant if the for- 
mation is instantaneous, while the core-Helium-burning 
stars may dominate the total NIR emission for a more 
constant star formation scenario. Therefore, the behavior 
described above is probably due to the distinct evolution 
of high-mctallicity AGB stars (see Willson 2000 and ref- 
erences therein) and core-Helium-burning stars (Mowlavi 
et al. 1998) compared with the relatively well-defined se- 
quence in their evolutionary properties established for sub- 
solar metallicities. Within the error intervals derived, the 
metallicity is mainly degenerate with the age, specially in 
those sets including [/-band data. This is probably due to 
the reduction in the age-extinction degeneracy thanks to 
the information provided by the [/-band data about the 
rest-frame UV. 

4.2.5. Stellar Mass 

The comparison between Figures 6b and 7b shows that 
the mean uncertainties in the iT-band mass-to-light ratio 
(or stellar mass) of intermediate-redshift galaxies are much 
lower than those derived for the nearby sample. First, it 
is important to note that in these figures we represent 
absolute errors. For a Solar-abundant 12Gyr-old nearby 



galaxy formed instantaneously the if-band mass-to-light 
ratio is ^1.3Mq/Lk,q, while for a 5Gyr-old galaxy at 
z=0.7 is ~O.4M /Lx,0. Therefore, the relative uncertain- 
ties, assuming the average absolute uncertainties given in 
Table 3 for AC rl =0.07mag, would be about 30 and 20 per 
cent, respectively for the nearby and intermediate-redshift 
galaxies. Although this still implies a significant improve- 
ment in the -FT-band mass-to-light ratio determination, it 
is also noticeable that the K filter now traces the rest- 
frame J-band luminosity, which is more affected by the 
misunderstanding about the actual IMF (see Bell & de 
Jong 2001). Finally, the J-band luminosity is also more 
sensitive to small differences between the assumed expo- 
nential star formation and the galaxy actual star formation 
history than the rest-frame -ftf-band data. 

As we pointed out in Section 4.2.2 the bias in the age de- 
termination toward lower age values also leads to a strong 
bias in the -FT-band mass-to-light ratio of intermediate- 
redshift galaxies due to the systemic decrease in the rest- 
frame J-band luminosity per Solar mass with the age of 
stellar population when the age is older than ~1 Gyr. 

4.3. High-redshift galaxies 

4.3.1. Formation Timescale 

With regard to the formation timescale, Figure 8a shows 
that the bias toward lower timescale values observed at 
intermediate redshift is even more pronounced at high- 
redshift. This bias is a natural consequence of the diffi- 
culty of deriving/predicting the long-term star formation 
history of a galaxy when it is still very young. This is also 
evidenced by the fact that the mean timescale uncertainty 
increases systematically with redshift for the same observ- 
ing errors and band-pass combinations. In Table 3 we also 
show that in many cases (BVRI vs. U+BVRI; SDSS vs. 
SDSS+2MASS) the mean timescale uncertainties increase 
when a larger number of observing bands is used, with a 
progressive reduction in this bias. As in the intermediate- 
redshift case the dominant degeneracy involving the galaxy 
formation timescale occurs with the age of the stellar pop- 
ulation. 

4.3.2. Age 

The large formation timescale uncertainty described 
above and the existence of a strong age-timescale degener- 
acy, specially at ages older than 100 Myr, lead to very large 
age uncertainties, even larger than those derived for the 
intermediate-redshift galaxies. The age-timescale degener- 
acy at ages younger than 100 Myr is significantly smaller 
because at these young ages a change in the formation 
timescale, which ranges between 200Myr-6Gyr (see Ta- 
ble 1), does not affect to the UV-optical-NIR colors of 
the stellar population. In other words, the degeneracy in 
timescale within this age range is complete and no corre- 
lation between the age uncertainty and any other uncer- 
tainty is expected. In this case the main degeneracies are 
the age-extinction and the age-metallicity ones. 

Moreover, the age-timescale degeneracy in combination 
with the bias in formation timescale described above are 
also responsible for the strong bias in age observed in Fig- 
ure 8a at ages older than ^50 Myr. The fact that the 
UVIJK set provides a better age and timescale determi- 
nation than the U+BVRI and U+BVRI+K sets demon- 
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strates the importance of obtaining JHK data in order 
to derive the properties of high-redshift galaxies. This is 
due to the fact that the JHK filters now cover the rest- 
frame V, R, and z' optical bands, where the changes due 
to galaxy evolution are more noteworthy and the infor- 
mation content about the galaxy properties is larger. In 
particular, the JHK filters would provide information si- 
multaneously about the presence of main-sequence stars 
(via the rest-frame V and R bands), core-Helium-burning 
stars (via the rest-frame V, R, and z' bands), and AGB 
stars (via the rest-frame z'-band; £>0.5Gyr). 

4.3.3. Dust Extinction 

Because of the extensive coverage of the UV range of the 
spectrum, the study of high-redshift galaxies using optical- 
NIR colors leads to very small dust-extinction uncertain- 
ties. In this sense, the dust-extinction uncertainties given 
in Table 3 at this redshift assuming an observing error of 
0.10 mag are in the range E(B - V r )=0.03-0.07mag. The 
dust extinction within the interval of uncertainty is mainly 
degenerate with the age of the stellar population. 

4.3.4. Metallicity 

Figure 8 shows that the metallicity uncertainty for the 
high-redshift sample does not show the strong metallic- 
ity dependence found in nearby and intermediate-redshift 
samples. Only when JHK NIR data are available the un- 
certainties at very high metallicities become significantly 
smaller than those derived for the low metallicity galaxies. 
As we commented in Section 4.2.4 for the intermediate- 
redshift case, this is probably due to the distinct signa- 
ture of high-metallicity core-Helium-burning stars (eg. in 
the number ratio of blue-to-red supcrgiants; Mowlavi et 
al. 1998) within the age range considered. In the case of 
a SSP galaxy, these stars dominate the rest-frame VRz' 
( JHK at 2=1.4) emission for ages younger than 0.4 Gyr, 
while the emission at shorter wavelengths comes mainly 
from main-sequence stars (see Chariot & Bruzual 1991). 
It is important to note that the core-Helium-burning stars 
may dominate the emission in the R and z' bands up to 
ages of 5 Gyr for larger formation timescales. Within the 
uncertainty intervals obtained, the metallicity is mainly 
degenerate with the age of the stellar population. 

4.3.5. Stellar Mass 

The if-band mass-to-light ratio uncertainties derived 
here are very small compared with those obtained from 
the nearby and intermediate-redshift samples, with values 
ranging between 0.01 and 0.06 M Q /Lk,q- If we adopt a K- 
band mass-to-light of O.27M q /Lk,0, which corresponds to 
the value expected for a 3 Gyr-old galaxy with Solar metal- 
licity, the relative uncertainty would range between 5% 
and 20%, depending of the band-pass combination avail- 
able. Figure 8b shows that there is also a strong depen- 
dence of the mass-to-light ratio uncertainty with the value 
of the mass-to-light ratio itself. In particular, a clear min- 
imum in its uncertainty is observed at ages older than 
8 Myr, which is probably associated with the evolution of 
the massive stars off the main sequence toward the red su- 
pergiant phase. During this part of the evolution a sudden 
change in the rest-frame z' luminosity and optical colors 
of a SSP is produced, which could explain why the uncer- 
tainty is particularly small around this age value. 



4.4. Effects of the Model Uncertainties 

In this section we analyze the results obtained when the 
optical- NIR colors of a sample of galaxies generated using 
the GISSEL99 models are compared with the predictions 
of the PEGASE evolutionary synthesis models. We have 
restricted this comparison to the nearby sample and the 
range of properties specified in Table 1. The results of this 
comparison are shown in Figure 9. 

4.4.1. Formation Timescale 

Figure 9a shows that the same bias toward lower values 
of the formation timescale that we noted for the inter- 
mediate and high-redshift samples is also present in this 
case (see Sections 4.2.1 & 4.3.1). The main reason for the 
existence of this bias is the small change in the optical- 
NIR colors of the stellar population with the timescale 
when the age t«r. Therefore, in order to compensate 
both the observing errors and the differences in the color 
predictions between the GISSEL99 and PEGASE models, 
the best-fitting solution has to be found at lower values 
of the timescale where the intrinsic dispersion of the col- 
ors is larger. The existence of this strong bias also leads 
to very small timescale uncertainties compared with those 
obtained using the GISSEL99 models. Within the uncer- 
tainty intervals derived, the dominant degeneracy involv- 
ing the galaxy formation timescale is the age-timescale de- 
generacy. 

4.4.2. Age 

With regard to the age determination, the uncertain- 
ties derived are very similar for the BVRI and U+BVRI 
sets. However, for those band-pass combinations includ- 
ing NIR data the ages derived are strongly biased toward 
younger ages. The reason for this bias, which also leads to 
significantly smaller age uncertainties, is the offset in the 
(J - H) and (H - K) colors between the GISSEL99 and 
the PEGASE model predictions (see Figure 1) due to the 
differences in the modelling of the upper RGB and AGB 
phases. In particular, Figure lb shows that the PEGASE 
models are ~0.07mag redder in (J — H) and ^0.04 mag 
redder in (H - K) than the GISSEL99 models within the 
age range 4-12 Gyr. Therefore, in order to compensate 
for this difference in color, the best-fitting solution usu- 
ally leads to younger ages, which within this age range 
imply bluer colors. Because the differences in the colors 
between the two models only occur in the NIR, the optical 
colors predicted by the PEGASE models at these younger 
ages should be bluer than those of the sample. There- 
fore, in order to compensate for this effect, the age bias 
described above has to be accompanied by strong biases in 
dust extinction and/or metallicity that would lead to red- 
der optical colors. Within the error intervals derived the 
total degeneracy is dominated by the age-timescale and 
age-extinction degeneracies. 

4.4.3. Dust Extinction 

As we commented above (see also Figure 9a) there is a 
strong bias in dust extinction estimates toward higher ex- 
tinction values when NIR data are used. This bias, along 
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with the metallicity bias described in Section AAA, re- 
sults in a global reddening of the optical colors but a small 
change in the NIR colors of the galaxies in the sample. On 
the other hand, at very high extinction values the uncer- 
tainties are also biased by the upper limit in E(B — V) im- 
posed during the data-models comparison procedure (see 
Table 1). The mean uncertainties derived, both in age, 
dust extinction, and metallicity are very similar to those 
obtained by using the GISSEL99 models. 

4.4.4. Metallicity 

The distribution of the uncertainty in metallicity shown 
in Figure 9a indicates that a strong bias toward higher 
metallicity values is present when NIR data are available. 
As we commented in Section 4.4.2, this bias is probably 
related with the age bias and the differences in the NIR 
colors predicted by the two sets of models. As in the case 
of GISELL 99 models, the comparison with the PEGASE 
models leads to a clear dependence of the uncertainty with 
the metallicity value itself, with smaller uncertainties at 
very high metallicities (see Section 4.1.4). 

4.4.5. Stellar Mass 

The results shown in Figure 9b with regard to the K- 
band mass-to-light ratios mainly reflect the biases in the 
galaxy property determination, with the stellar masses de- 
rived systematically smaller than the input values. This 
is due (1) to the bias toward younger ages described in 
Section 4.4.2 and (2) to the higher if-band luminosity per 
unit mass of the PEGASE models compared with the GIS- 
SEL99 models (see Figure lb). Because of the stronger 
bias in age, the mean uncertaintines in the -fT-band mass- 
to-light ratio are smaller than those obtained using the 
GISSEL99 models (see Section 4.1.5). Finally, Figure 9b 
shows that the mass-to-light ratio uncertainty becomes 
higher at older ages and lower timescale values. This be- 
havior, which is also present in the case of the GISSEL99 
models (see Section 4.1.5), is due to the progressive in- 
crease in the rate of change of the if -band mass-to-light 
ratio with log t (see Figure 1) accompanied by a small de- 
pendence of the age uncertainty (in log t scale) with the 
value of the age itself. 

5. CONCLUSIONS 

In this study we have analyzed the dependence of the un- 
certainties and degeneracies in the galaxy properties upon 
different parameters: (1) the combination of bands avail- 
able, (2) the observing errors, and (3) the galaxy proper- 
ties themselves (including redshift). 

Here we summarize our main results and point out some 
directions for the optimization of galaxy evolution studies 
using broad-band photometry data. We describe sepa- 
rately the nearby, intermediate, and high-rcdshift cases. 

Nearby galaxies: In order to determine the star for- 
mation history, age, and dust extinction of nearby galaxies 
with relatively small uncertainties the use of U -band data 
is fundamental. The availability of if-band data also al- 
lows a reduction in the uncertainty in the age and metal- 
licity of the stellar population, but the use of additional 
J and f/-band data is largely redundant. The use of the 
i-T-band data is unfortunately limited by the existence of 



large uncertainties in the modeling of the if-band lumi- 
nosities and NIR colors of stellar populations. The most 
significant reduction in the age and if-band mass-to-light 
ratio uncertainty is achieved when UV data are used. The 
poor treatment of the post-AGB and "extreme" HB phases 
by the existing evolutionary synthesis models introduce, 
however, an additional uncertainty during the data-model 
comparison, which is particularly important in the case of 
very old stellar populations. For the same number of ob- 
serving bands, the availability of wider wavelength base- 
lines results in lower uncertainties. Both the formation 
timescale and i-T-band mass-to- light ratio uncertainties are 
larger when the corresponding values for these properties 
are larger. On the other hand, the metallicity uncertainty 
decreases with the value of the metallicity itself due to 
the distinct photometric evolution of high-metallicity core 
Helium burning stars. 

A complete description of the physical reasons behind 
these conclusions and of the degeneracies responsible for 
the uncertainties described above are given in Section 4.1 
(see also Section 4.4). 

Intermediate-redshift galaxies: The star formation 
history of intermediate-redshift galaxies can be derived 
with worse precision than in nearby galaxies because their 
stellar populations are younger. The age uncertainty is 
smaller than in the nearby-galaxies case and shows a 
strong bias toward younger ages. A significant reduction of 
this bias and of the mean uncertainties is achieved when 
NIR data are used, especially if all three J, H, and K- 
band data are available. The dust-extinction uncertainty 
is larger for larger values of the dust extinction itself. The 
use of U -band data provides an important reduction of this 
dependence and of the mean dust-extinction uncertainty. 
If U -band data are available the use of additional UV data 
do not provide much more information about the galaxy 
properties. The use of NIR data (J, H , and if-band data) 
significantly reduces the uncertainty in the metallicity of 
the galaxy. The absolute and relative uncertainties in the 
galaxy if-band mass-to-light ratio are smaller than those 
derived for nearby galaxies. However, the fact that the K 
filter now covers the rest-frame J-band leads to a larger 
uncertainty associated with the IMF and with the param- 
eterization of the galaxy star formation history and, con- 
sequently, to a larger stellar mass uncertainty. For a more 
detailed description see Section 4.2. 

High-redshift galaxies: As expected, the bias and 
mean uncertainty in the determination of the timescale 
for the galaxy formation are even larger in this case that 
in the nearby or intermediate-redshift galaxies. The age 
of the stellar population is derived with a large uncer- 
tainty, only reduced when JHK data are available. The 
dust-extinction in these galaxies can be derived to a very 
high accuracy even when only optical data are available. 
The use of JHK data is fundamental in order to improve 
both the age and metallicity determinations. A complete 
description of the uncertainties and degeneracies between 
these properties is given in Section 4.3. 

Some of the conclusions drawn above can also be found 
through the literature expressed in a qualitative way. How- 
ever, this work constitutes the first systematic and quan- 
titative study on the optimization of broad-band photom- 
etry for studies on the evolution of galaxies. It is impor- 
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tant to note that the application of these results to future 
galaxy surveys can help to reduce the uncertainty in the 
derivation of the galaxy physical properties, sometimes by 
weakening a particular degeneracy but most of the time by 
decreasing the intervals over which this degeneracy takes 
place. In Table 3 we have summarized the mean uncer- 
tainties in the galaxy stellar population properties derived 
in this paper considering different redshifts, sets of observ- 
ables, and observing errors. 

Our results are directly applicable to spectrophotomet- 
ric surveys like the SDSS and surveys looking for emission- 
line galaxies at fixed redshifts (Martin, Lotz & Ferguson 
2000; Moorwood, van der Werf, Cuby & Oliva 2000; Iwa- 
muro et al. 2000; Pascual et al. 2001; Zamorano et al., in 
preparation). However, in the case of the blind- rcdshift 
surveys a comparison between our results and those from 
previous studies on the optimization of the photometric- 
rcdshifts technique (Kodama et al. 1999; Bolzonella et al. 
2000; Mobasher & Mazzei 2000; Wolf et al. 2001) is still 
needed. 

6. FUTURE APPLICATIONS 

The results summarized above demonstrate that the de- 
sign of galaxy evolution studies based only on qualitative, 
intuitive ideas may lead (in some cases avoidably) to large 
uncertainties. 



Because of this we intend to apply this work to the de- 
sign of future projects on galaxy evolution estimating opti- 
mal sets of observables and required signal-to-noise ratios. 
Although in this work we have only considered broad-band 
filters, this procedure is easily generalizablc to combina- 
tions of broad, medium, and narrow-band filters from the 
far-UV to the near-infrared. In addition, the combined 
use of the procedure here described with state-of-the-art 
radiative transfer and dust models (Popescu et al. 2000) 
will allow us to extend this range up to sub-millimeter 
wavelengths. 

Beyond the results shown in this paper, we can also de- 
rive, upon request, the uncertainties and degeneracies in 
the galaxy properties for a given combination of filters, 
observing errors, and galaxy redshift. 
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Fig. 1. — Comparison of the predictions of the Bruzual & Chariot (in preparation; thick-lines), PEGASE (version 2.0; see Fioc & Rocca- 
Volmcrangc 1997; thick grey lines), and Starburst 99 (Lcitherer et al. 1999; thin-lines) models for a Solar mctallicity galaxy with instantaneous 
(SSP; panel a) and continuous star formation (panel b). 
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Fig. 2. — a, b & c) Top panels: Distribution of solutions associated with the 200 Monte Carlo simulations in the times cale- age space for 
a nearby galaxy with a formation timescale of 4.5 Gyr, an age of 5 Gyr, a dust extinction E(B — V)=0.08 mag, and Solar mctallicity. Middle 
panels: Distribution of solutions in the age-metallicity space for the same galaxy model. Bottom panels: The same for the age- extinction 
space. Observing errors are 0.10, 0.07, and 0.03 mag, respectively for the a, b, and c figures. The input properties of the galaxy are marked 
with a star symbol. This comparison was performed for the GALEX+SDSS+2MASS color set. The size of each point is proportional to the 
value of the maximum likelihood estimator, C, for the corresponding Monte Carlo simulation. Degeneracies between the galaxy properties 
are evidenced by correlations in the distribution of points. 
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Fig. 3. — Uncertainties in the derived properties as function of the properties themselves for a subsamplc of 500 nearby galaxies and the 
U+BVRI+JHK set of observables. In Panel a the mean differences between the derived and the input properties computed in intervals of 
0.5 Gyr in formation timescalc, 0.025 mag in color excess, and 0.05 dex in age are represented by a light-grey line. Mean ±1(7 values for the 
derived properties are delimited by two dark-grey lines. The region defined by the mean ± 1 <j lines is also represented in Panel b. Grey 
shaded areas represent the regions not covered by our comparison procedure. 




Table 1 

Range in galaxy properties from which synthetic galaxy colors of the sample were generated and range in the 

model parameters for the data-models comparison. 



Sample Models 



E(B - Thenar 0.00-1. 00 mag 0.00-1.05 mag 

r 0.2-6 Gyr 0.01-10 Gyr 

Age (z=0) 4-12 Gyr 1-15 Gyr 

(2=0.7) 1-5 Gyr 0.5-8 Gyr 

(2:=1.4) 3Myr-3Gyr lMyr-6Gyr 

Metallicity (Z) 0.2, 0.4, 1.0, 2.5, 5xZ 0.2, 0.4, 1.0, 2.5, 5xZ 
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Fig. 4. — a) Frequency histograms for the components of the PCA1 vectors for the sample of high-rcdshift galaxies with AC„=0.07mag 
for three different sets of observables (SDSS, SDSS+2MASS, and GALEX+SDSS+2MASS). This figure suggests that the age-metallicity and 
age-extinction degeneracies are comparable for the SDSS set, while the age-extinction and the age-timescale degeneracies are competing in 
the case of the SDSS+2MASS and GALEX+SDSS+2MASS sets, b) Cosine of the angle between the PCA1 vector and the age-timescale 
(top), age-extinction (middle), and agc-metallicity (bottom) planes, as measure of the degeneracy between the stellar populations properties 
(see Section 4), as funcion of the age for high-rcdshift galaxies with an uncertainty AC n =0.07mag and the GALEX+SDSS+2MASS set. The 
sign of the cosine indicates if the degeneracy is in the sense that an increase in both properties can lead to the same observational properties 
(positive) or the value for one magnitude has to be decreased while the other is increased (negative). From this figure is clear the dependence 
of the dominant degeneracy with the age of the stellar population, going from the age-metallicity and age-extinction degeneracies at ages 
below 10 Myr, to only the age-extinction at ages between 10 Myr and 300 Myr, and the age-timescale at ages between 300 Myr and 3 Gyr (see 
Section 4.3.2). 
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Table 2 

Definition of the different combinations of filters analyzed 



Bands 


# of bands 


Alias 


B,V,Rc,Ic 


4 


BVRI 


U,B,V,Rc,Ic 


5 


U+BVRI 


U,B,V,R c ,Ic,K 


6 


U+BVRI+K 


U,B,V,R c ,Ic,J,H,K 


8 


U+BVRI+JHK 


B,V,R c ,Ic,J,H,K 


7 


BVRI+JHK 


B,V,Rc,Ic,K 


5 


BVRI+K 


U,V,I C ,J,K 


5 


UVIJK 


u', g', r', i', z' 


5 


SDSS 


u',g',r',i',z',J,H,K 


8 


SDSS+2MASS 


FUV.NUV.u', g', r', %' , z', J, H, K 


10 


GALEX+SDSS+2MASS 
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Fig. 5. — Uncertainties in the derived properties of nearby galaxies for the BVRI, U+BVRI, and U+BVRI+K sets (panel a) and the SDSS, 
SDSS+2MASS, and GALEX+SDSS+2MASS sets (panel b) assuming observing errors of 0.03 mag. 
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Fig. 6. — Uncertainties in the derived if -band mass-to-light ratios for the BVRI+K, U+BVRI+K, and U+BVRI+JHK sets (panel a) and 
the BVRI+K, SDSS+2MASS, and GALEX+SDSS+2MASS sets (panel b) assuming observing errors of 0.03 mag. 
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Fig. 8. — a) Uncertainties in the derived properties of high-rcdshift galaxies for the BVRI, U+BVRI, and UBVRI+K sets assuming 
observing errors of 0.03 mag. b) _ft"-band mass-to-light ratio uncertainties for the BVRI+K, U+BVRI+K, and U+BVRI+JHK sets assuming 
observing errors of 0.10 mag. 
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Fig. 9. — Mean timescale, age, dust extinction, metallicity (panel a), and A"-band mass-to-light ratio (panel b) uncertainties derived using 
the PEGASE evolutionary synthesis models. In panel a we show the results for the BVRI, U+BVRI, and UBVRI+K sets, and in panel b 
those for the BVRI+K, U+BVRI+K, and U+BVRI+JHK sets. In both cases observing errors of 0.03 mag have been adopted. 
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Table 3 

Mean Ict uncertainties in the derived properties 
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Property 


Unit 




2=0 






2=0.7 
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mag 


0.07 


0.14 


0.16 


0.02 


0.04 


0.05 


0.02 


0.03 


0.04 




r 


Gyr 


2.06 


2.65 


2.89 


2.36 


3.07 


3.34 


3.20 


3.41 


3.39 




log t 


dex 


0.25 


0.36 


0.39 


0.09 


0.16 


0.21 


0.12 


0.21 


0.26 




log Z 


dex 


0.08 


0.19 


0.25 


0.10 


0.16 


0.21 


0.08 


0.17 


0.23 




(M/L)jr 


M Q /Lx,0 


0.24 


0.35 


0.37 


0.02 


0.04 


0.05 


0.01 


0.01 


0.02 


GALEX+SDSS+ 




mag 


0.04 


0.09 


0.12 


0.01 


0.02 


0.02 


0.01 


0.02 


0.03 


+2MASS 


T 


Gyr 


1.51 


2.26 


2.60 


2.01 


2.93 


3.26 


2.92 


3.63 


3.67 




log t 


dex 


0.10 


0.19 


0.24 


0.08 


0.13 


0.17 


0.09 


0.16 


0.21 




log Z 


dex 


0.05 


0.16 


0.22 


0.11 


0.16 


0.20 


0.03 


0.08 


0.14 




(M/L) K 


Mq/Lx, 


0.10 


0.23 


0.28 


0.02 


0.03 


0.04 


0.01 


0.01 


0.01 



